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ABSTRACT 

We use a volume-, magnitude-limited sample of nearby galaxies to investigate the 
effect of the environment on the Hi scaling relations. We confirm that the Hi-to-stellar 
mass ratio anti correlates with stellar mass, stellar mass surface density and NUV — r 
colour across the whole range of parameters covered by our sample (10^ ;SM»;S10^^ M©, 
7.5 ^/z*;£9.5 Mq kpc~^, 2 i>NUV — ri:6 mag). These scaling relations are also followed 
by galaxies in the Virgo cluster, although they are significantly offset towards lower 
gas content. Interestingly, the difference between field and cluster galaxies gradually 
decreases moving towards massive, bulge-dominated systems. By comparing our data 
with the predictions of chemo-spectrophotometric models of galaxy evolution, we show 
that starvation alone cannot explain the low gas content of Virgo spirals and that only 
ram-pressure stripping is able to reproduce our findings. Finally, motivated by previous 
studies, we investigate the use of a plane obtained from the relations between the Hi- 
to-stellar mass ratio, stellar mass surface density and NUV — r colour as a proxy 
for the Hi deficiency parameter. We show that the distance from the 'Hi gas fraction 
plane' can be used as an alternative estimate for the Hi deficiency, but only if carefully 
calibrated on pre-defined samples of 'unperturbed' systems. 

Key words: galaxies: evolution-galaxies: fundamental parameters-galaxies: clus- 
ters: individual: Virgo-ultraviolet: galaxies- radio lines:galaxies 



1 INTRODUCTION 

Understanding the role played by the cold gas component in 
galaxy evolution is one of the main challenges for extragalac- 
tic studies. The cold atomic hydrogen (Hi) is the reservoir 
out of which all stars are eventually formed. Therefore a de- 
tailed knowledge of its relation to other galaxy properties is 
of primary importance in order to build a coherent picture 
of galaxy evolution. 

Although the direct detection of Hi emission from 
individual galaxies is still technically l imited to the 
nearb y universe (e.g., IVerheiien et~aLll2007l : [ Catinella et al.l 
|2008| ). in the last decades several studies have high- 
lighted how the Hi content (usually estimated as the 
Hi mass-to-luminosity ratio) varies with morpho logy, lu- 
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to be the use of 21 cm data to investigate environmen- 
tal effects. Several studies have shown that Hl-deficien10 



galaxies are mainly /only pre sent in v ery high-density en- 
vironments llGiovanelli fc Hay ncs 1985j_ Cavatte et al.lll990l : 
ISolanes et al. I I2OOII : iBoselli fc Gavazzil 120091 '). providing 
strong constraints on the possible environmental mecha- 
nisms responsible for the quenching of the star fo rmation 
in clusters of galaxies (see lBoselli fc Gavazzil [2OO6I and ref- 
erences therein). 

Unfortunately, a great part of what we know about Hi in 
galaxies has been achieved by studying relatively small sam- 
ples, with not always well defined selection criteria, often fo- 
cused on specific morphological classes (e.g., only late-type 
galaxies) and with small overlap with multiwavelength sur- 
veys, necessary to investigate how the cold gas correlates 
with other baryonic components. Moreover, it is very diffi- 
cult for simulations to provide accurate estimates for impor- 
tant observables like the Hi deficiency parameter (since it is 
based on a detailed morphological classification), limiting 
the comparison between observations and theory. In other 
words. Hi astronomy is in a situation similar to the one 
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^ The Hl-deficiency {Defjji) is defined as the difference, in loga- 



rithmic units, between the expected Hi mass for an isolated galaxy 
with the same morphologica l type and optical diameter of the tar- 
get and the observed value llHavnes &: Giovanelli|[T984h . 
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of optical and ultraviolet astronomy before the Sloan Digi- 
tal Sky Survey and the GALEX mission. The main scaling 
relations (e.g., colour-magnitude, size-luminosity, etc.) were 
already known, but a detailed quantification of their proper- 
ties, crucial for a comparison with models, was still missing. 

Luckily, the situation is rapidly changing for 21 cm stud- 
ies. The advent of large Hi surveys, such as the Arecibo 
Legacy Fast ALFA Survey (|Giovanelli et al.l 12005'). which 
eventually will provide Hi masses for a few tens of thousands 
of galaxies over areas of sky with large multiwavelength cov- 
erage, is gradually allowing detailed statistical analysis of 
Hi properties in the local universe. Very promising results 
in this direction have been rece ntly obtained by the GALEX 
Arecibo SDSS Survey (GASS: [Catinella et allboiOl ). a tar- 
geted survey of a volume-l imited sample o f ~100 massive 
galaxies (M. >W^° Mq). ICatinella etall (|2010l ) used the 
first GASS data release to quantify the main scaling rela- 
tions linking the Hl-to-stellar mass ratio to stellar mass, stel- 
lar mass surface density and colour. They also suggested the 
existence of a 'gas fraction plane' linking Hi gas fractior0 to 
stellar mass surface density and NUV — r colour and pro- 
posed to use this plane to isolate interesting outliers that 
might be in the process of accreting or loosing a significant 
fraction of their gas content. In this context, the plane might 
be considered as an alternative way to define Hi deficiency 
when accurate morphological classification is not available. 
However, this hypothesis has never been tested for a sample 
for which both Hi deficiency and the Hi gas fraction plane 
can be estimated independently. 

In this paper, we use a volume-, magnitude-limited sam- 
ple of '^SOO galaxies to extend the study of the Hi scaling 
relations to a larger dynamical range in stellar mass and 
stellar mass surface density and to different environments 
(i.e., from isolated systems to the center of the Virgo clus- 
ter). Our main goals are a) to accurately quantify the effect 
of the cluster environment on the Hi scaling relations, b) to 
investigate the validity of the Hi plane as a proxy for Hi de- 
ficiency and c) to highlight the power of Hi scaling relations 
to constrain models of galaxy formation and evolution. 



2 THE SAMPLE 

The analysis presented in this paper is based o n the Herschel 
Reference Survey (HRS, i Boselli et all [2OI0I '). Briefly, this 
consists of a volume-limited sample (i.e., 15< d <25 Mpc) 
including late-type g alaxies (Sa and later) with 2MASS 
l|Skrutskie et al.ll2006l ) K-band magnitude Kstot < 12 mag 
and early-type galaxies (SOa and earlier) with Kstot < 8.7 
mag. Additional selection criteria are high galactic latitude 
(b > -|-55°) and lo w Galactic extinction (A_b < 0.2 mag, 
ISchlegel et al.|[T998l ). to minimize Galactic cirrus contami- 
nation. The total sample consists of 322 galaxies (260 late- 
and 62 early -type galaxies) 0. As extensively discussed in 
iBoselli et al.l (2010), this sample is not only representative 
of the local universe but, spanning different density regimes 



In this paper we will refer to the M (HI)/Mt ratio as gas frac- 
tion. 

^ HRS228/LCRS B123647.4-052325 has been removed from the 
original sample since the redshift reported in NED was not correct 
and the galaxy lies well beyond the limits of the HRS, z ~0.04. 



(i.e., from isolated systems to the center of the Virgo clus- 
ter), it is also ideal fo r environmental studies (see also 
ICortese fc Hughedl2009l l . 

Atomic hydrogen masses have been estimated from 
Hi 21 cm line emission data (mai nly single-dish) , avail- 
able from the literature (e.g.. ISpringob et all 20051: 
iGiovaneUi et al.ll2007l : iKent et al.ll2008l : iGavazzi et al.ll2003l 
and the NASA/IPAC Extragalactic Database, NED). In to- 
tal, 305 out of 322 galaxies in the HRS (~95%) have been 
observed at 21 cm, with 265 detections and 40 non detec- 
tions. Hi masses have been computed via 

^ = 2.356 X 10% I"' W 
Mq Jy kms-i VMpc/ ^ ' 

where Shi is the integrat ed Hi line flux-densi ty and d is the 
distance. As discussed in iBoseUi et all (|2010l ). we fixed the 
distances for galaxies belonging to the Virgo cluster (i.e., 
23 Mpc for th e Virgo B clo ud and 17 Mpc for all the other 
clouds; Gavaz zi et al.|[l999h . while for the rest of the sample 
distances have been estimated from their recessional veloci- 
ties assuming a Hubble constant Hq =70 km s"'^ Mpc~^. 
In case of non detections, upper limits have been deter- 
mined assuming a 5a signal with 300 km s~^ velocity width. 
We estimate the Hi deficien cy parameter {Defni) following 
iHavnes fc Giovanellil (|l984l '). The expected Hi mass for each 
galaxy is determined via 

\og{M{HI),^p) = am + hni x log (^J^) - 21og(fe) (2) 

where h — Hq/IQQ km s~^ Mpc~^, D25 is the optical isopho- 
tal diameter measured at 25 mag arcsec^ in B-band and aui 
and bni are two coefficients that vary with morphological 
type. These coefficients have been calc ulated for the f ollow- 
ing types (se e also Table 3 in BoscUi &: GavazzillioO^ 'l: SOa 
and earlier llHavnes &: Giovanellil 11984 ) , Sa-Sab, Sb, Sbc, 
Sc |Solanes et al! 19961 ) and Scd to Irr (jBoselli fc Gavazzi 
120091 ) . Morphological classifications fo r our sample are take n 



from the Vi rgo Cluster Cat alogue (|Binggeli et al.l 119851 ) . 
NED and .Bo selli et al.1 (|2010l ). The Hi deficiency is then 



DefHi = log(A/(i//)e 



log(M(iJJ)) 



(3) 



We note that the estimate of the expected Hi mass for early- 
type galaxies is extremely uncertain. As we will show later, 
this uncertainty can significantly affect the quantification of 
the Hi scaling relations. In the following, we will consider 
as 'Hl-deficient' galaxies those objects with Defni > 0.5 
(i.e., galaxies with 70% less hydrogen than isolated sys- 
tems with the same diameter and morphological type). 
The average Hi deficiency of galaxies with Defni < 0.5 
is Z)e/H/=0.06±0.28, consistent with the typical uncer- 
tainty in the estimate of Defni for isolated fi eld galaxies 
(jHavnes fc Giovanelli|[l984l : ISolanes et al.lll996l ). 

Ultraviolet and optical broad-band photometry have 
been obtained from Galaxy Evolution Explorer (GALEX, 
iMartin et"al] l2005h and Sloan D igital Sky Survey DR7 
(SDSS-DR7. lAbazaiian et al.ll2009ll databases, respectively. 
GALEX near-ultraviolet (NUV; A=2316 A: AA=1069 A) 
images have been mainly obtained as part of two on- 
going GALEX Guest Investigator programs (GI06-12 , P.I. 
L. Cortese and GI06-01, GUViCS. iBoselU et aLllioTH) . Ad- 
ditional frames have been obtained from the GALEX GR6 
public release. All frames have been reduced using the cur- 
rent version of the GALEX pipehne (ops-v7). The GALEX 
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Figure 1. The Hi mass fraction of the HRS plotted as a function of stellar mass, NUV — r colour, concentration index and stellar 
mass surface density. Black circles and green triangles represent detections and non-detections, respectively. GASS DRl detections and 
non-detections are indicated with magenta dots and arrows, respectively. The average scaling relations obtained from ALFALFA stacking 
of the GASS parent sample are indicated in cyan. 



NUV and SDSS g,r,i photometry for the HRS was de- 
termined as follows. The SDSS images were registered to 
the GALE X frames and convol ved to the NUV resolu- 
tion (5.3", iMorrissev et al.l |2007| ). Isophotal ellipses were 
then fit to each image, keeping fixed the center, ellip- 
ticity and position angle (generally determined in the i- 
band). Asymptotic magnitudes have been determined from 
the growt h curve obtained follo wing the technique de- 
scribed bv lGil de Paz et al.l(l2007l) and corrected for Galac- 
tic extinction assuming a " iCardeUi et al] l| 19891 ) extinction 
law with A(V)/E( B - V) =3.1: i.e., A{X)/E{B - V)=8.2 
l|Wvder et al. 1'2007). 3.793, 2.751 and 2.086 for NUV, g, r 
and i, respectively. Stellar masses M* are determined from 
i-band luminosities Li using the g — i colour-dependent stel- 
lar mass- to-light ra t io rela tion from lZibetti et all (|2009| ). as- 
summg a lChabrieil (|2003l ) initial mass function (IMF): 

log(M,/LO = -0.963 + 1.032 * {g - i) (4) 



The final sample used for the following analysis includes 
those HRS galaxies for which Hi, NUV and SDSS obser- 
vations are currently available: 241 galaxies (~75% of the 
whole HRS, namely 200 late- and 41 early- type galaxies). 



3 THE Hi SCALING RELATIONS OF THE HRS 

In order to quantify how the Hi gas fraction varies as a func- 
tion of integrated galaxy properties, we plot in Fig[T]the Hl- 
to-stellar mass ratio as a function of stellar mass (upper-left 
panel), observed NUV — r colour (upper-right), concentra- 
tion index in r-band (C3i(r), defined as the ratio between 
the radii containing 75% and 25% of the total r-band light) 
and stellar mass surface density (i.e., Mt/(27rii|o,i)i where 
-Rso.i is the radius containing 50% of the total i-band light 
i.e., the effective radius). For compar ison, we also plot the 
data from the GASS data release 1 (|Catinella et al.ll2010l . 
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Figure 2. The Hi mass fraction as a function of stellar mass, NUV — r colour and stellar mass surface density for Hl-normal galaxies 
(upper panel), galaxies outside the Virgo cluster (middle panel) and galaxies belonging to one of the Virgo cluster clouds (bottom panel). 
Symbols are as in Fig. [T] 
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Figure 3. The average Hi mass fraction (i.e., < log{M(H I) /M*) >) scaling relations for different samples. Hl-normal, galaxies belonging 
to or outside Virgo are indicated by blue, red and green lines respectively. For comparison, GASS DRl is shown in magenta. 



magenta symbols) and from the stacki ng of ALFALFA ob- 
servations for the GASS parent sample (|Fabello et al.|[201ll . 
cyan line). Over the range of stellar masses covered by the 
HRS, the Hi gas fraction strongly anti-correlates with M, , 
NUV — r (a proxy for specific star formation rate, SSFR) 
and fj,, , while a very weak non-linear trend is observed with 
the concentration index. The tightest correlation is with the 
observed NUV — r colour (Pearson correlation coefficient 
r — —0.89, dispersion along the y-axis a =0.42 dexjf], while 
the scatter gradually increases for the stellar surface density 
(r = -0.74, a =0.61 dex) and mass (r = -0.72, a =0.63 
dex). 

Ou r findings extend very nicely the results of the GASS 
survey (|Gatinella et al.ll2010l : iFabello et al.ll201ll ) to lower 
stellar mass and surface densities, highlighting how the star 
formation rate, stellar mass and galaxy structure are tightly 
linked to the Hi content of the local galaxy population. How- 
ever, it is important to notice that the scaling relations for 
the HRS are slightly offset towards lower gas fractions when 
compared to the results of the GASS survey. This difference 
is easily understood when we consider that, by construc- 
tion, nearly half of the galaxies in the HRS belong to the 
Virgo cluster and, as a whole, this sample might be biased 
towards gas-poor objects. To test this and to characterize 
the Hi scaling relations in different environments, we have 
divided our sample in three different subsets: (a) galaxies 
with Defni < 0.5 (i.e., Hl-normal, 135 galaxies), (b) galax- 
ies outside the Virgo cluster (110 galaxies) and (c) galaxies 
belonging to one of the Vir go cluster clou ds (Virgo A, B, 
N, E and S as defined bv .Gavazzi et al . 1999, 131 galaxies). 
We note that, while (b) and (c) are complementary by con- 
struction, (a) may include both Virgo and field galaxies. The 
scaling relations for the three samples are plotted in Fig. [21 
As expected, although for all samples the Hi gas fraction 
decreases with stellar mass, colour and stellar mass surface 



^ These parameters are obtained assuming the non-detections to 
their upper-limits. 



density, galaxies in different environments show significantly 
different Hi content. Virgo galaxies have, on average, a lower 
Hl-to-stellar mass ratio than Hl-normal galaxies across the 
whole range of stellar masses and stellar mass surface den- 
sities covered by the HRS. In fact, the dispersion in the two 
scaling relations drops from ~0.62 dex for the whole sample 
to ~0.37 dex for Hl-normal galaxies only. 

In order to quantify the difference between cluster 
and field, in Fig. [3] we show the average trends (i.e., < 
log{M{HI)/M,) >) for each subsample. The averages are 
measured by assuming the non detections to their upper- 
limit and by correcting for the fact that we do not have 
Hi and ultraviolet data for all HRS galaxies. In practice, 
we weighted each galaxy by the completeness of the sample 
in its bin of K-band luminosity. The average scaling rela- 
tions and the number of galaxies contributing to each bin 
are given in Table [T] 

The difference between Hl-normal and cluster galaxies 
varies by ~0.3 dex as a function of stellar mass (i.e., from 
~0.55 to 0.85 dex), while it increases from ~0.25 to 1 dex 
when moving from disk- to bulge-dominated systems (i.e., 
from /i, ~ 7.4 to 8.7 Mq kpc^^jf]- 

Less strong is the en- 
vironmental dependence of the gas fraction vs. NUV — r 
relation. Blue-sequence galaxies {NUV — r <3.5 mag) have 
the same Hi gas fraction regardless of the environment, con- 
sistent with a scenario in which actively star-forming cluster 
galaxies have just started their infall into the cluster center 
(|BoseUi fc Gavazzill200d ICortese et al.ll2008al ). On the con- 
trary, red Virgo galaxies are significantly gas poorer than 
Hl-normal systems. It is important to note that this differ- 
ence might be, at least partially, due to an extinction effect. 



^ We remind the reader that, by construction, the HRS does 
not include intermediate-, low-mass ear ly-type galaxies, typically 
found in high-density environments llHughes &: Cortesj |2009| : 
ICortese fc Hughesll2009l) . Thus, for low stellar masses, low stellar 
mass surface densities and red colours, we might be underesti- 
mating the difference between cluster and field galaxies. 
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While a significant fraction of red cluster galaxies are likely 
to be really passive systems, field red objects should be more 
heavily affected by internal dust attenuation. Therefore bins 
of observed NUV — r colour might not directly correspond 
to bins of SSFR. We plan to investigate this issue in a future 
paper, where Herschel data will be used to properly correct 
for internal dust absorption. 

Interestingly, we find some differences also when 
comparing Hl-normal galaxies to objects outside the 
Virgo clusters. Since Hl-deficier it galaxies are mainly /only 
found in clusters of galaxies ( Havnes fc Giovanellj Il984l : 



Table 1. The average scaling relations for the three samples dis- 
cussed in Sec. 3. 



ICortese et al ] |2008bl : [Boselli fc G avazzi 200§), it IS common 
practice to assume that galaxies in the field, in pairs and 
loose groups are not Hl-deficient (and viceversa). Our anal- 
ysis shows that a one-to-one correlation between Hl-normal 
galaxies and systems in low-/intermediate-density environ- 
ments apparently breaks at high stellar masses (A/» >10^'''^ 
Mq), stellar mass surface densities (/x, > 8.5 Mq kpc~^) and 
red colours {NUV — r >4.5 mag), where we find Hl-deficient 
galaxies also outside the cluster environment. It is important 
to note that a significant fraction of these gas-poor galaxies 
(i.e., 5 out of 7 and 6 out of 10 for the log(M(i?/)/M.) vs. 
log(M«) and log(M(i//)/M«) vs. log(/i*) relations, respec- 
tively) are early-type systems (SOa or earlier) for which the 
Hi deficiency calibration is notoriously uncertain. Whether 
the typical red, massive, bulge-dominated galaxy contains 
a significant cold gas reservoir or not is still matter of de- 
bate. These ranges of stellar masses, colour and stellar sur- 
face densities have not been adequately sampled by previous 
Hi investigations and it is not completely surprising that we 
find a difference between Hl-normal and field galaxies. 

Whatever the origin of the difference between Hl- 
normal and 'field' galaxies, the break observed in the 
Hi scaling relations for galaxies outside the Virgo cluster 
is very similar to the transition observed in the colour- 
stellar mass and colour-stellar ma s s surface densi t y dia - 
gra ms (e.g.. iKauffniann et al.ll20o3 ). iBothwell et~al] ()2009l l 
and ICatinella et al.l I 2010l ) alreadv reported the presence 
of a transition in the Hi gas fraction vs. stellar mass and 
stellar mass surfa ce density relations, respectively (see also 
iKannappanI |2003 ). Our analysis confirms both results at 
once independently, suggesting that the typical transitions 
observed in the colour-mass diagram might be related to the 
transition in Hi gas fraction showed in Figs. [2] and |3] Un- 
fortunately, above the transition our statistics becomes too 
small and we will have to wait for larger Hi surveys to char- 
acterize more in detail this break in the Hi scaling relations. 

Finally, it is important to point out that, above the 
transition in stellar mass/surface density, the difference in 
Hi gas fraction between cluster and field galaxies tends to 
disappear. This provides an indirect support to a scenario in 
which environmental effects are today mainly active on the 
population of low-/intermediate-mass disk galaxies, while 
the properties of high-mass/bulge-dominated systems are 
less affected by the local density. 



4 COMPARISON WITH MODELS 

The results presented in § [3] have shown that, at fixed 
stellar mass and surface densities, cluster galaxies have 
significantly lower Hl-content than galaxies in low den- 
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a: We note that t hese are averages of lo g(M {H I) / Mt) and not 
M{HI)/M, (as in lCatinella et al.|[2oTol) . We preferred this ap- 
proach because the distribution of Hi gas fraction is closer to a 
log-normal than to a Gaussian. 

sity environments. On one side, these results confirm that 
Hl-deficient ga laxies are mainly present in high density 
environments l|Havnes fc Giovanellil 1 1984 ). On the other, 
they nicely complement previous analysis showing that 
cluster galaxies have in general a lower star formation 
rate than their field counterpart (iBoselli fc Gavazzil |2006| : 
iGavazzi et"ai]|2002l :[G omez et al.ll2003l )." suggesting that the 
gas los s is behind the quen c hing of the star fo r matio n ac- 
tivity (|Gavazzi et all l2006l : ICortese fc Hughes! |2009| ). Al- 
though this scenario is commonly accepted, there is still 
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Figure 4. Model predictions for a starvation (upper panels) and ram-pressure stripping (lower panels) event. In case of starvation we 
considered ages between and 9.5 Gyr, while for ram-pressure we restricted our analysis to events occurred between and 5 Gyr ago. 
In both sets of panels, the unperturbed model is indicated by the filled circles. In the bottom panel, the two lines indicate models for 
different ram-pressure efficiencies: eo =0.4 (light curve) and 1.2 (dark curve) Mq kpc~^ yr~^- Data are as in Fig. [3] We note that these 
are not evolutionary tracks, but illustrate where a galaxy would lie today depending on the age of the interaction. 



debate about the environmental mechanism (s) responsible 
for such diff erences between low- a nd high-density environ- 
ments (e.g.. IWeinmann et al.|[201ll ). In particular, it is not 
clear whether the Hi is stripped directly from the disk via 
a strong interaction with t he dense intra-cluste r medium 
(i.e., ram-pressure stripping: iGunn fc Gott"l972'). or if the 
gas is just removed from the extended gaseous halo sur- 
roun ding the galaxy, pr e venting furth e r in fall (i.e., starva- 
tion; [LMSonet^lli^. [BoseiirelaD (|2008i ) have recently 
used multi-zone chemical and spectrophotometric models to 
constrain the evolutionary history of dwarf galaxies in the 
Virgo cluster and to discriminate between the ram-pressure 
and starvation scenarios. They showed that the properties 
of most of the dwarfs dominating the faint end of the Virgo 
luminosity function are consistent with a scenario in which 
these galaxies were initially star-forming systems, accreted 



by the cluster and stripped of their gas by one or more ram 
pressure stripping events. 

Here, we us e the same approach adopted by 
iBoselli et ahl l|2008D to carry out a similar test for intermedi- 
ate and massive galaxies in the Virgo cl uster. The details of 
the models here adopted can be found in iBoselU et al. I (|2008l . 
Appendix B). Briefly, the evolution of galaxies is traced us- 
in g the multi-zone chemical and spectrophotometric model 
of lBoissier fc PrantzosI (|200d ). updated with an em pirically 
determined star formation law (|Boissier et al.|[2003l ) relating 
the sta r formation rate to t he total gas surface densities. Fol- 
lowing iBoselli et al. I ('2008'), we modified the original models 
to simulate the effects induced by the interaction with the 
cluster environment. In the starvation scenario we simply 
stopped the infall of pristine gas in the model. We investigate 
interactions that started between and 9.5 Gyr ago. The 
ram pressure stripping event is simulated by assuming a gas- 
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loss rate inversely proportional to the potential of the galaxy, 
with an efficiency dependi ng on the density of the intra- 
cluster medium (taken from lVollmer et al]l200ll ). We assume 
two different values for the peak efficiency of ram-pressure 
stripping (eo =0.4 and 1.2 M© kpc^^ yr~^) and an age of the 
interaction between and 5 Gyr ago. Models have been gen- 
erated assuming a spin parameter A =0.05 and four different 
rotational velocities (Vc =100, 130, 170, 220 km s"^), in or- 
der to cover the whole range of stellar masses spanned by the 
HRS. In order to convert the output of the model into the 
observables here considered, we proce eded as follows. Ste l- 

jKroupaet al.l(|l993h . 



lar masses have been converted from a 



used in the model, to alChabrieil (120031') I MF by adding 0.06 



dex (|Bell et al.ir2003l : ICallazzi et al.ll200i 'l. SteUar mass sur 



face densities have been determined using the effective ra- 
dius in i-band, as done for our data. Hi masses are obtained 
from the total gas mass by removing the contribution of he- 
lium and heavy elements (i.e., a factor 1.36) and assuming a 
molecular-to-a tomic hydrogen gas ra tio of 0.38, independent 
of stellar mass l|Saintonge et al.ll201ll ). NUV ~r colours have 
been 'reddened' by assuming an average internal extinction 
A(NUV) =1 mag. This value is quite arbitrary and will 
likely depend on the evolutionary history of each galaxy. 
However, this assumption does not affect our main conclu- 
sions. Finally, we note that the goal of this exercise is to es- 
tablish which scenario is more consistent with out data, not 
to determine the best model fitting our observations. There 
are many free parameters and assumptions in the modeling 
that could be tweaked in order to better match our observa- 
tions, but looking for an exact fit does not seem meaningful 
given the simple description adopted here. 

In Fig. [3] (top panels) we compare the predictions of 
the starvation model with our data. The filled circles indi- 
cate the predictions for the unperturbed models. The models 
seem to fairly reproduce the NUV — r and vs. gas fraction 
scaling relations observed for Hl-normal and field galaxies. 
However, a significant difference between data and models 
is observed in the A4{HI)/AI, vs. Al, relation: for very high 
stellar masses (rotation velocities), the unperturbed model 
predicts gas fractions ^^0.5 dex higher than the average ob- 
served value. This likely reflects the fact that the model is 
only valid for disk galaxies and does not include any bulge 
component. It is thus not surprising that at high masses, 
where the fraction of bulge-dominated galaxies increases, 
the predictions we obtain are not representative of the whole 
population (although still well within the range of gas frac- 
tions covered by the observations). This is also confirmed by 
the fact that the model covers just the range of stellar mass 
surface densities typical of disk galaxies. 

By comparing the predictions of the starvation model 
with our data, it clearly emerges that starvation is only able 
to mildly affect the Hi content, even for very old events. 
The Hi gas fraction decreases by only ~0.25 dex, i.e. within 
the intrinsic dispersion of the scaling relation for Hl-normal 
galaxies. This decrease is even lower in the M{HI) / vs. 
NUV — r plot, where galaxies affected by starvation lie on 
the main relation observed for unperturbed galaxies. It is 
thus impossible to reproduce the difference between cluster 
and field galaxies by just stopping the infall of gas onto the 
disk. 

Much more promising are the results obtained for the 
ram-pressure stripping model (Fig. |4l bottom panels). For 
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Figure 5. Left panel: The Hi gas fraction plane for Hl-normal 
{Defni <0.5) HRS galaxies. The Hi gas fraction plane obtained 
for GASS galaxies is indicated by the magenta line. Right panel: 
The correlation between Hi deficiency and distance (along the 
y-axis) from the plane. The solid line indicates the best linear 
fit to the Hi detections. In both panels, Hl-normal, Hl-deficient 
galaxies and non-detections are indicated with black filled, red 
empty circles and green triangles, respectively, and the dotted 
line shows the 1:1 anti correlation. 



all the three scaling relations here observed, ram-pressure 
stripping is able to fairly reproduce the difference observed 
between cluster and field galaxies. A maximum efficiency of 
0.4 Me kpc"^ yr"^ (light curve) is sufficient to reproduce 
the average decrease in Hi content for Virgo galaxies, while 
a higher efficiency (1.2 Mq kpc~^ yr~^, dark curve) is nec- 
essary to explain some of the most Hl-deficient systems in 
our sample. As expected, the M{H I)/M, vs. NUV — r is 
the relation least affected by the environment. Once the gas 
is removed, the star formation is quickly reduced and the 
perturbed galaxy will not shift considerably from the main 
relation. In conclusion, only the direct removal of gas from 
the star-forming disk is able to reproduce the significant dif- 
ference between the Hi scaling relations of cluster and field 
galaxies. 



5 THE Hi GAS FRACTION PLANE AND THE 
Hi deficiency PARAMETER 

In the previous sections we have confirmed the existence of 
a tight relation between Hi gas fraction, stellar mass surface 
density and NUV — r colour. Recent works have used these 
three quantities to define a 'gas fraction plane' to be used 
in order to deter mine gas fractions f or large samples lacking 
Hi observations (jZhang et all |2009| . using optical colours) 
and to isolate galaxies in the process of accreting or loosing 
a signific ant amount of gas yZJatincUa ct al. 2010). As dis- 
cussed bv lZhang et al.l (|2009l ). the existence of this plane can 
be seen as a direct consequence of the Kennicutt-Schmidt re- 
lation. In order to test how good the gas fraction plane is 
as a proxy for Hi deficiency, we decided to compare the two 
approaches for our sample. We thus fitted a plane to the two- 
dimensional relation between Hi mass fraction, stellar mass 
surface density and NUV — r colour for Hl-normal galaxies 
only (Defni <0.5), following the methodology outlined in 
I Bernard! et al.l l|2003h . i.e., by minimizing the residuals from 
the plane. The best fit to the plane is 



log 



M{HI) \ _ 
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-0.33(iV(/V -r)~ 0.40 log(/i,) -K3.73 (5) 
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Figure 6. Left panel: The Hi gas fraction plane for HRS galaxies 
outside the Virgo cluster. The Hi gas fraction plane obtained for 
Hl-normal HRS galaxies (Fig. 4) is indicated by the solid line. 
Right panel: The correlation between Hi deficiency and distance 
(along the y-axis) from the plane. Symbols are as in Fig. [5] 

and has a dispersion of only ~0.27 dex, comparable to 
the typical uncertainty in the Hi deficiency parameter. The 
plane so obtained is illustrated in Fig. [5] (left panel) and 
is fairly co nsis tent with the one determined for the GASS 
sample bv lCatinella et al.l (|201()l . magenta line). Given the 
different datasets and selection criteria adopted by the two 
surveys, this agreement is quite remarkable and indicates 
that both Hl-normal HRS and GASS samples provide a fair 
representation of the Hi properties of massive galaxies in the 
local universe. As expected, Hi-deficient galaxies (empty cir- 
cles) are significantly offset from the plane, having a lower 
Hl-content than what expected for their colour and stel- 
lar surface density. The distance from the plane along the 
y-axis (i.e., log{MHi / Mt)obs - log {Mr i /M,) plane) strongly 
anti correlates with the Hi deficiency (Fig. [5l right panel). 
The best linear fit for the Hi detections is 

Defni = (-1.09 ± 0.04) x Aplane + (0.14 ± 0.02) (6) 

with ~0.25 dex dispersion. Again, the scatter is similar to 
the intrinsic scatter in the definition of the Hi deficiency, 
suggesting that the Hi plane is a tool as good as the Hi de- 
ficiency to isolate extremely Hl-rich or Hl-poor galaxies, in 
particular when accurate morphological classification is not 
available. 

Although the results obtained here are very promising, 
we want to conclude this section with some notes of cau- 
tion in the definition, and use, of the Hi gas fraction plane. 
The definition of a plane for Hl-normal galaxies is justified 
by the fact that, for this sample, we observe a linear corre- 
lation between log(M(_H'/)/M,) and both NUV - r colour 
and log(/i,). As discussed in the previous section, these lin- 
ear relations are not valid for typical 'field' galaxies across 
the whole range of parameter here investigated, but we find 
a break for massive, bulge-dominated, red galaxies. There- 
fore, when these objects are included, the concept of a plane 
is no longer justified. In order to show how important is the 
sample selection, we plot in Fig. [6] the plane obtained for 
galajcies outside the Virgo cluster and compare it with what 
previously obtained for Hl-normal galaxies (solid line). The 
two planes are significantly different, in particular at low 
Hi gas fractions where there is an offset up to a factor ~4 
in the estimate of the Hi mass. Incidentally, this shows once 
more how 'field' galaxies is not the same as 'Hl-normal' sys- 
tems. Moreover, this new plane is almost useless to isolate 



Hl-deficient galaxies. Hl-poor systems are no longer strong 
outliers and the relation between Hi deficiency and distance 
from the plane is too scattered (cr ~0.4 dex) to be useful. 
We note that the good agreement between our plane for Hl- 
normal galaxies and the one obtained using GASS detections 
is due to the fact that both datasets sample the stellar mass 
surface density regime below the break in Hi gas fraction. 
Notice that, if the GASS detection limit was much deeper 
than 1.5% gas fraction (which roughly corresponds to the 
limit separating Hl-normal and Hl-deficient galax;ies at high 
stellar masses, see Fig.[2|, then this survey would include de- 
tections above the /i, break, where the definition of a plane 
is no longer justified. 

In conclusion, the Hi gas fraction plane is a promising 
tool to estimate gas fractions and isolate outliers, but only if 
calibrated on well defined samples of 'normal' unperturbed 
galaxies. 



6 CONCLUSIONS 

In this paper we have shown the relations between Hl- 
to-stellar mass ratio, galaxy structure and colour for the 
HRS, a volume-, magnitude-limited sample of nearby galax- 
ies. Hi gas fractions anti correlate with stellar mass, stellar 
mass surface den sity and NUV — r col our, confirming previ- 
ous results (e.g.. ICatinella et al.ll2O10h . We show that these 
trends are observed in both low- and high-density environ- 
ments, with the only difference that the scaling relations 
for cluster galaxies are shifted towards lower gas fractions. 
Thus, the cluster environment seems to have just a sec- 
ondary (although important) role on galaxy evolution, with 
the efficiency of gas consumption and star formation mainly 
drive n by the intrins i c 'size' (e.g., total mass) of the system 
(G ayazzi et al.l Il996l : ISoselh et all l200ll : iKauffmann et al.l 
I2OO3I ). Our findings nicely complement the scaling rela- 
tions between SSFR and galaxy pr operties revealed by 
large optical and ultrav iolet surveys (|Blanton et al.l l2005l : 
ISchiminovich et al.ll2007t ). Naturally, because the Hi is the 
fuel for star formation, it is likely that the link between 
SSFR and integrated properties is driven by the scaling re- 
lations presented here. This is also supported by the presence 
of a break in the M{HI)/Mt vs. ^* and M{HI)/Mt vs. M* 
relations, very close to the typical value characterizing the 
transition between the blue cloud and red sequence. 

In order to shed light on the environmental mechanism 
responsible for the lower Hi content of cluster galaxies, we 
compared our results with th e predictions of the chemo- 
spectrophotometric model by iBoissier fc PrantzosI l|2000h . 
We find that simply halting the infall of pristine gas from 
the halo, mimicking the effect of starvation, is not sufficient 
to reproduce our observations. Only the stripping of gas di- 
rectly from the star-forming disk (i.e., ram-pressure) is able 
to match our data. This naturally explains why Hl-deficient 
galaxies are mainly/only found in very high-density environ- 
ments, where ram-pressure stripping is efficient. However, it 
is important to note that our results do not imply that star- 
vation is not playing any role on galaxy evolution. They just 
show that, if galaxies suffer of starvation, their star forma- 
tion activity and gas content can be changed, but the shape 
of the main Hi scaling relations will not be significantly al- 
tered. We remind the reader that these results are only valid 
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for disk galaxies. At high masses, where bulge-dominated 
systems dominate, not only our model cannot be used but, 
most importantly, the difference between field and cluster 
starts to disappear, suggesting that the origin of the lower 
Hi content of these objects cannot only be due to the cluster 
environment. 

One of the interesting outcomes of our analysis is that, 
at high stellar masses, 'field' galaxies are not the same as 
'Hl-normal' objects. If confirmed, this clear difference might 
have important implications in our understanding of envi- 
ronmental effects in the local universe. In particular, one of 
the main issues for environmental studies is to reconcile the 
evidence that the effects of the environment start at de nsi- 
ties typical of poor groups (|BIanton &: Moustakas|[2009l and 
references therein) with the fact that Hl-deficient objects are 
only observed in clusters of galaxies. Our analysis suggests 
that, from a statistical point of view, Hl-deficient galaxies 
are clearly segregated in high-density environments only for 
M*JS10^° M© and /i.SlO* M© kpc"^ or, in other words, 
that intermediate-/low-mass gas-poor quiescent disks are 
only found in clusters of galaxies. Thus, the results obtained 
from Hi data turn out to be fully consistent with what shown 
by optical studies: i.e., the intermediate-/low-mass part of 
t he red sequence i s only populated in clusters of galaxie s 
jBaldrv et al.ll2006l : iHaines et ai]|2008l : iGavazzi et al1l2010f ) . 

Finally, we have confirmed that the Hi gas fraction 
plane defined bv lCatinella et al.l (i2010) can be used as an al- 
ternative to the Hi deficiency parameter. However, we have 
also illustrated how important is the sample selection in the 
definition of the Hi gas fraction plane and that, as for the 
Hi deficiency, its use is mainly valid for galaxies below the 
threshold in /i,. 

Our results clearly highlight the importance of Hi scal- 
ing relations for environmental studies and the necessity to 
extend this approach to larger samples and to a wider range 
of environments in order to understand in detail the role 
played by nurture in galaxy evolution. 
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